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ABSTRACT 

We present threshold enhanced QCD corrections to rapidity distributions of di-leptons in the Drell- 
Yan process and of Higgs particles in both gluon fusion and bottom quark annihilation processes 
using Sudakov resummed cross sections. We have used renormalisation group invariance and the 
mass factorisation theorem that these hard scattering cross sections satisfy as well as Sudakov 
resummation of QCD amplitudes. We find that these higher order threshold QCD corrections 
stabilise the theoretical predictions under scale variations. 



Perturbative Quantum Chromodynamics (pQCD) provides a framework to successfully com- 
pute various observables in the collisions of hadrons at high energies. Recent theoretical advances 
in the computations of higher order QCD radiative corrections have lead to precise results for sev- 
eral important observables. Because of this progress, we can now predict these observables with 
unprecedented accuracy for physics studies at the Tevatron collider in Fermilab as well as at the 
upcoming Large Hadron Collider (LHC) in CERN [1]. 

The Drell-Yan (DY) production of di-leptons [2] has been one of the most important probes of 
the structure of hadrons. It is also one of the dominant production processes at hadron colliders. 
At the LHC, it will serve as a luminosity monitor which is very important to precisely calibrate the 
machine for searches for physics beyond the Standard Model (SM). In DY production, a pair of 
leptons is produced through the decay of virtual photons, Z and W bosons that result from the colli- 
sions of incoming partons (quarks and gluons) in the hadrons. At hadron colliders, the DY process 
provides precise measurements of various standard model parameters. Rapidity distributions of 
Z bosons [3] and charge asymmetries of leptons coming from W boson decays [4] can probe the 
structure of the hadrons and possible excess events in di-lepton invariant mass distributions can 
point to physics beyond the standard model such as i?-parity violating supersymmetric models and 
models with Z', or with contact interactions [5]. Both DO and CDF collaborations [6] at the Fermi- 
lab Tevatron made precise measurements of Z and W production cross sections and asymmetries 
which not only allowed for stringent tests of the standard model but also play an important role in 
the Higgs search at future colliders. These measurements are also possible at the LHC due to the 
large cross sections for the DY process. 

The other process which is equally important is Higgs boson production at these colliders 
because it will establish the Standard Model as well as look for beyond the SM Higgs [7, 8]. The 
Higgs boson, which is responsible for the electroweak symmetry breaking in the Standard Model, 
is yet to be discovered. The search for this particle has been going on at the Fermilab Tevatron and 
is one of the most important tasks for the CERN LHC. The LEP experiments in the past provided 
vital information on the possible mass range of this particle [9]. The lower bound on the mass is 
1 14.4 GeV/c 2 and an upper bound is less than 219 GeV/c 2 at 95% CL [10]. At the LHC, Higgs 
bosons will be predominantly produced through the gluon fusion process due to the large flux of 
gluons in the protons at these energies. They can be detected through the rare two photon decay 
mode which has less QCD background than other signals. 

In pQCD, the total cross sections for the DY production of di-leptons and Higgs boson pro- 
duction are known upto next-to-next-to-leading order (NNLO) level [11-27]. However, due to the 
complexity involved with the top quark loops, the Higgs production cross sections are only known 
in the large top quark mass limit beyond the next-to-leading order (NLO). In addition to these fixed 
order results, the resummation programs for the threshold corrections to the total cross sections for 
DY and Higgs production have also been very successful [28,29](see also [30]). See [22,31] for 
next- to-next-to-leading logarithmic (NNLL) resummation results. Due to several important QCD 
results at the three loop level that are recently available [32-38], the resummation upto N 3 LL has 
also become a reality [39-42]. The fixed order partial soft-plus-virtual N 3 LO corrections [39,43] 
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to the DY and Higgs productions show the reliability of the perturbation theory and the stability 
against the scale variations. The fixed order results as well as the resummed results reveal very 
interesting structures in the perturbative QCD results (see, [42,44^-7]). 

Infra-red safe observables, such as hadronic cross sections, are computed using the QCD im- 
proved parton model. Due to the factorisation property that certain hard scattering cross sections 
satisfy, they can be expressed interms of finite partonic cross sections convoluted with the parton 
distributions functions (PDFs). The partonic cross sections are calculated in QCD using standard 
perturbation theory in powers of the strong coupling constant g s that becomes small at high en- 
ergies. The ultraviolet singularities that arise beyond leading order are often removed in the MS 
renormalisation scheme at a renormalisation scale /jr. The collinear singularities that result due 
to the presence of light mass partons are mass factorised into the bare parton densities in the MS 
scheme at a factorisation scale /jf . Hence the fixed order perturbative results are often sensitive 
to these scales /jr and hf- For example see [48] for a study of scale variations of the Higgs cross 
section in NLO. However results that are known to sufficiently high order in the strong coupling 
constant are often less sensitive to these scales because the observables are renormalisation group 
invariant. That is they would be strictly independent of choice of these scales if the entire pertur- 
bative expansion were known. 

In addition to the scale uncertainties the fixed order computations suffer from the presence of 
various large logarithms which arise in some kinematical regions. These regions are often impor- 
tant from the experimental point of view. The large largarithms spoil the standard perturbative 
approach. The alternate approach is to resumm these logarithms in a closed form. Resumming a 
class of large logarithms supplemented with fixed order results can usually cover the entire kine- 
matic region of the phase space. In this paper, we will mainly concentrate on a class of logarithms 
that arise in the threshold regions. These threshold corrections are further enhanced when the 
fluxes of the incoming partons become large in those regions. In the case of Higgs production 
through gluon fusion, the gluon flux at small partonic energies becomes large enhancing the role 
of threshold corrections. Here we examine the effects of soft gluons that originate in the threshold 
region of the phase space when we consider the xf and rapidity distributions in DY production 
and Higgs production through both gluon fusion and bottom quark annihilation. Here the large 
logarithms are generated when the gluons that are emitted from the incoming/outgoing partons 
become soft. 

In [42], we found that soft distribution functions of Drell-Yan and Higgs production cross sec- 
tions in perturbative QCD are maximally non-abelian. That is, we found that the soft distribution 
function for Higgs production can be obtained entirely from the DY process by a simple multipli- 
cation of the colour factor Ca/Cf- In the article [43], using the soft distribution functions extracted 
from DY, and the form factor of the Yukawa coupling of Higgs to bottom quarks, we predicted 
the soft-plus-virtual (sv) part of the Higgs production through bottom quark annihilation beyond 
NNLO with the same accuracy that the DY process and the gluon fusion to Higgs process are 
known [39,42]. We extended [44] this approach to entirely different processes such as Higgs de- 
cay to bottom quarks and hadroproduction in e + e annihilation. The approach that we followed 
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in [42-44] is closely related to that of the standard threshold resummation and hence we could de- 
termine [42] the threshold exponents D\ upto three loop level for DY and Higgs production using 
our resummed soft distribution functions and B\ for both deep inelastic scattering and Higgs decay 
and hadroproduction. In this paper we extend this approach to include differential cross sections 
such as x F and rapidity distributions of the di-lepton pair in DY production and of Higgs bosons in 
Higg production processes. 

In the following we systematically formulate a framework to resum the dominant soft gluon 
contributions to these differential cross sections. We perform the resummation in the Zi(i =1,2) 
space of the kinematic variables, which are the appropriate scaling variables that enter the differ- 
ential partonic cross sections. The threshold region corresponds to z, — > 1 and in this region all 
the partonic cross sections are symmetric in z\ <-> zi- We have used renormalisation group (RG) 
invariance, mass factorisation and Sudakov resummation of QCD amplitudes as guiding principles 
to perform the resummation in this region. Using the resummed results in z\ space we predict 
the soft-plus-virtual parts (also called threshold corrections) of the dominant partonic differential 
cross sections beyond N 2 LO. We also study the numerical effect of our predictions on both the 
xp and rapidity distributions of di-leptons and Higgs bosons. The analytical results are presented 
in the Appendices for both DY and Higgs production 1 through gluon fusion. For an early refer- 
ence where the resummation for DY differential distributions at rapidity Y = (or xp = 0) was 
considered consult [49]. 

The differential cross section can be expressed as: 
do 1 

^ = <^orn(4Aq 2 )W I (x° l ,xiq 2 ), I = q,b,8 , (D 

with the normalisation W^ om {x Q v x^q 2 ) = 8(1 — Xj)8(l — x®). The x® (i = 1,2) are related to the 
kinematical variables q 2 and x. Here q is the momentum of the di-lepton pair in the DY process 
and of the Higgs boson in the Higgs production. The variable x can be the xp or rapidity of the 
di-lepton pair or of the Higgs boson. For di-lepton production, I = q and G 7 = d(5 q (x,q 2 ,x) /dq 2 
with q 2 the invariant mass of the final state di-lepton pair i.e, q 2 = Mf +l _ . For Higgs production 
through gluon fusion, I = g and o 7 — <5 8 (%,q 2 ,x) and for Higgs production through bottom quark 
annihilation I = b and G 7 = o h (x,q 2 ,x) with q 2 = mjj where m H is the mass of the Higgs boson. 
The variable x = q 2 /S with S = {p\ + P2) 2 the center of mass energy squared where are the 
momenta of incoming hadrons Pi (i = 1,2). In the QCD improved parton model, the function 
W^x^x^q 2 ) can be expressed in terms of the PDFs appropriately convoluted with perturbatively 
calculable partonic differential cross sections denoted by A ! d ab , where the subscript "d" stands for 
"differential", as follows 

W 1 {x\,X Q 2 ,q 2 ) = £ / dxi / dx 2 Xa b (x h X 2 ,l4) 

lr rhe results for the Higgs production via bottom quark annihilation are not presented here but can be obtained 
from the authors on request. 
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x f dzi [ dz2 5(x? -XiZl) 8(X2-X 2 Z2) &d,ab(zi, Z2, ■ (2) 

Jo Jo 

Here, [xr is the renormalisation scale and hf the factorisation scale. We consider the differential 
cross sections for two kinematic variables namely 

i{pi-pi)-q A v !, / P2-q\ 

x = x F = , and x = Y = -log . (3) 

S 2 \p\-qj 

For the xf (x = x F ) distribution, the x® variables satisfy 

Xf=X° 1 -X2, X = x\x2 , (4) 

while for the rapidity Y (x = Y) distribution, we have 

Here, the function Hj lh {x\,X2^ F ) is the product of PDFs f a (xi,/ii F ) and fb(x2,Mp) renormalised at 
the factorisation scale fip. That is, 

^ b (x h x 2 ,fi 2 F ) = fPixuflp) f^ 2 (x 2 ,/4), 

X* b (xi,X2,i4) = x 1 f^{x h /4)x 2 f b 2 {x2,i4), (6) 
with xi (i= 1 , 2) the momentum fractions of the partons in the incoming hadrons. 

We first study the contributions coming from the soft gluons. The infra-red safe contributions 
from the soft gluons can be obtained by adding the soft part of the differential cross sections with 
the ultraviolet renormalised virtual contributions and performing mass factorisation using appro- 
priate counter terms. This combination is called the "soft-plus- virtual" (sv) part of the differential 
cross section. We call the remaining part the hard part of the differential cross section. Hence we 
write 

A d,ab(zhZ2,q 2 ,/4iVl) = A S(^1^2,^ 2 ,A'F,A'i)+S^A s ^ / (zi,Z2,^ 2 ,A'F,A'i), I = q,b,g. (7) 

The contributions coming from the hard parts A^(zi,z2,<? 2 ,A , F 5 A'i) of the differential cross sec- 
tions can be obtained by the standard procedure discussed in detail in [50, 51]. The soft-plus-virtual 
parts of the differential cross sections (A s ^ / (zi,z2,^ 2 ,A ,2 ?7A'f)) are found to be 

AS d,i(zhZ 2 ,q 2 ,M 2 R ,l4') = C exp ( x ¥ I d (q 2 ,^i 2 R ,fi 2 F ,zi,Z2,e) 



(8) 

£=0 



where ^ I d (q 2 ,/Li\,iLi 2 T ,zi,Z2,£) are finite distributions They are computed in 4 + £ dimensions and 
take the form 

V d (q 2 ,H 2 R A,zi,Z2,z) = ^ln(z / (4,//i// 2 ,e)) 2 + ln|F / (a„(2 2 , / / 2 ,e)|^5(l-z 1 )5(l-z 2 ) 

+2 ^>/(4,^ 2 ,// 2 ,zi,z 2 ,e)-clnr / /(4,/i 2 ,// 2 7 ,zi,e) 8(1 -z 2 ) 

-C \nT u {a s ,n 2 , ^F,Z2,e) 8(1 -zi), I = q,b,g. (9) 
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The symbol "C" means convolution. For example, C acting on the exponential of a function 
f(z\,zi) means the following expansion: 

Ce f(zi,z 2 ) = 8 (i_ Zl ) 8 (i_ Z2 ) + i. / ( Zl>Z2 ) + i. / ( Zl)Z2 ) (g)/ ( Zl>Z2 ) 

+^f{zi,z 2 ) ®/(zi,z 2 ) ®/(zi,z 2 ) + • • • • (10) 
In the rest of the paper the function f{z\ , zi) is a distribution of the kind 8(1— Zj) or £> ; - (z 7 -) , where 

-\n\l-Zj) 



i = 0,l,..-, and ;=1,2, (11) 



and the symbol ® means the "double" Mellin convolution. It convolutes with respect to the vari- 
ables z\ and Z2 separately. Since we are only interested in the sv part of the cross sections, we drop 
all the regular functions that result from various convolutions. F 7 (a^, <2 2 ,jU 2 , e) are the form factors 
that contribute to di-lepton (7 = q) (in DY) and Higgs (/ = g, b) production. In the form factors, we 
have Q 2 = —q 2 . The partonic cross sections depend on two scaling variables zi,z 2 . The functions 
4>J(a J ,<7 2 ,^ 2 ,zi,Z2,£) are called the soft distribution functions. The unrenormalised (bare) strong 
coupling constant d s is defined as 



8 s 



— 16*- (12) 

where g s is the strong coupling constant which is dimensionless in n = 4 + £ space time dimensions. 
The scale n comes from the dimensional regularisation which makes the bare coupling constant g s 
dimensionless in n dimensions. The bare coupling constant d s is related to renormalised one by 
the following relation: 

fu 2 ^ 

2\„ t.2\ / ^ 



Sea s = Z(fj R )a s (fj R ) , (13) 

where 5e = exp { | [y E — \n4-n] } is the spherical factor characteristic of n-dimensional regularisa- 
tion. The renormalisation constant Z(//|) relates the bare coupling constant d s to the renormalised 
one a s {^ R ). They are both expressed in terms of the perturbatively calculable coefficients |3 ; which 
are known up to four-loop level in terms of the colour factors of SU(N) gauge group: 

N 2 -l 1 
C A =N, Cf = -^, 7> = -. (14) 

Also we use rif for the number of active flavours. In the case of Higgs production, the number of 
active flavours is five because the top degree of freedom is integrated out in the large m top limit. 

The factors Z / (a. s ,/j|,/j 2 ,£) are the overall operator renormalisation constants. For the vector 
current Z q (d s ,iu Rl fj 2 1 e) = 1, but both the gluon operator [52] and the bottom quark coupling to 
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Higgs [53] get overall renormalisations. They satisfy the following RG equations: 




renormalised form factors F g and F respectively. 

The bare form factors F 1 {d s ,Q 2 ^ ,z) of both fermionic and gluonic operators satisfy the fol- 
lowing differential equation that follows from the gauge as well as the renormalisation group in- 
variances [54-57]. In dimensional regularisation, 



where the K contain all the poles in £ and the G collect the rest of the terms that are finite as 
£ becomes zero. The fact that the F 1 <2 2 ,,u 2 ,£) are renormalisation group invariant and the 
functions G 1 are finite implies that the K 1 terms can be expressed in terms of finite constants A 1 , 
the so-called cusp anomalous dimensions and the coefficients p\. The solution to the eqn.dToT) can 
be obtained as a series expansion in the bare coupling constant in dimensional regularisation. The 
formal solution up to four-loop level can be found in [34, 42] . 

The boundary conditions on the Sudakov differential equation, denoted by G/(e) (see eqn.(19) 
of [43]) can be found for both I = q and / = g in [35] to the required accuracy in £. We have 
extended this in [43] to the form factor corresponding to the Yukawa interaction of the Higgs 
boson to the bottom quarks. These constants G/(e) are expressed in terms of the functions B\ and 
fj. The B\ are known up to order cr s through the three-loop anomalous dimensions (or splitting 
functions) [32, 33] and are found to be flavour independent, that is Bf = i?f . The constants fj are 
analogous to the cusp anomalous dimensions A\ that enter the form factors with A? = A\ . It was 
first noticed in [27] that the single pole (in £) of the logarithm of the form factors upto two-loop 
level (a 2 ) can be predicted 2 due the presence of constants fj because these f 1 are found to be 
maximally non-abelian obeying the relation 



similar to the A\. This relation has been found to hold even at the three loop level [35]. With 
this information we can now predict all the poles of the form factors at every order in d s from 
these constants A 1 ,B l ,f , their anomalous dimensions, and the finite parts of the lower order (in a s ) 
contributions to the form factors. Interestingly, the single pole terms in the form factors contain 
the combinations [27,34,42,43] 




(16) 




?g 

i ' 



(17) 




2 A similar analysis of the structure of single pole terms of four-point amplitutes at the two-loop level can be found 



in [58,59]. 
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at order a' s . The terms — 28/, g i p\_i — 28/ ; ^_ 1 come from the ultraviolet divergences that are 
present in the loop integrals. These pole terms go away when the form factors undergo over- 
all operator UV renormalisation through the renormalisation constants Z 1 which satisfy the RG 
equations given in eq. (fT51) . 

The collinear singularities that arise due to massless partons are removed using the mass fac- 
torisation kernel T(zj,^,B) in the MS scheme (see eqn.®). We have suppressed the dependence 
on d s and in T. The factorisation kernel T(zj,/iip,E) satisfies the following renormalisation group 
equation: 

d , , 1 



UfUp £ 



(18) 



where the P(zj,/J F ) are the well-known DGLAP matrix- valued splitting functions which are known 
upto three-loop level [32,33]: 



P(zj,M 2 F ) = ^M)P 
(=1 



(19) 



The diagonal terms in the splitting functions P^ l \zj) have the following structure 



+ P, 



reg,// 



(zj) 



(20) 



where jj(zj) are regular when the argument approaches the kinematic limit (here Zy — » 1). The 
RG equations can be solved by expanding in powers of the strong coupling constant. For the soft- 
plus-virtual part of the differential cross sections, only the diagonal parts of the kernels contribute. 
We find the solutions contain only poles in e in the MS scheme: 



r(zy,A&e)=8(l-zy) + £4 



f4_ 



Zj,£) 



(21) 



The constants rW(zy,e) expanded in negative powers of £ up to four- loop level can be found 
in [42]. The Tu{d s ,jji 2 ,jjp,Zj,£) in eqn.© is the diagonal element of T(zj, /dp, e). 



The fact that A S J are finite in the limit £ 



implies that the soft distribution functions should 
have a pole structure in £ similar to that of F 1 and Tjj. To systematically study the soft distribution 
functions, we demand that they satisfy similar Sudakov type differential equations that the form 
factors F 1 satisfy (see eqn.dToT)): 



d 



tf,i?, z\ , Z2, e) 



q 2 jj 2 R 

&si ~ T i ~ T)ZljZ2>£ 

f4 r 



(22) 



where again the constants Kj contain all the singular terms in £ and the gJ are finite functions of 
£. Also the functions 4>J(a. s ,^ 2 ,^ 2 ,zi,Z2,e) satisfy the renormalisation group equations: 



A—^^d{d s ,q 2 ^ 2 ,zuZ2^)=Q. 



(23) 



7 



This renormalisation group invariance leads to the following equations 



d r'f \ A 



-rl, 



rijjfXd [**, ^ZUVL,* ) = -A K(^))5(l -zi)8(l -Z2 



d -^if „ q 2 /J 2 R 



ft R ft 



-rl . 



fii^Ga[a s ^ 1 %z h z 2 ^)=A\a s (4))8(l-zi)8(l-Z2). 



(24) 



If <&j(d s ,q ,/j ,zi,Z2,&) contains the correct poles to cancel the poles coming from F ,Z and Tu 
in order to make A s ^ l finite, then the constants A 1 have to satisfy 



1' = -^. 



(25) 



Using the above relation, the solution to the RG equation for G d [d s ,q 2 / iji 2 r ,ij? r / ij 2 ,z\,Z2-,t) is 
found to be 

,2 



\ a s,-T,-<2,Zl,Z2,£ 
fi R f 1 



Gd ( ^G4),-T,zi,Z2,e 



G d ( a s(q 2 )A>ZhZ2,e) 

-5(1 -zi)5(l -Z2) 2 ^ (a,(^ 2 )) . (26) 



7/4 

With these solutions, it is straightforward to solve the Sudakov differential equations yielding 

4>J(4,<?V, £1,22, e) = 4> ( /(a,,^ 2 (l-z 1 )(l-z 2 ), A / 2 ,£) 
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i=i 



, (g \l- Zl )(l- Z2 ) 
ft 1 



(i e) : 



4(l-zi)(l-z 2 ) 



where 



^' (,) (e)+Gi' ( °(8) 



$i' (0 (e)(27) 



(28) 



The above solutions for satisfy the fact that A^ v 7 are finite as £ — > (see eqn.©). The constants 
Kj (e) are determined by expanding in powers of the bare coupling constant d s as follows 



K l d [d s , p,zi,z 2 ,e ) = 5(1 -zi)5(l -z 2 ) 



i,(4V T ^zi'« ( E) 



(=i 



A' 



(29) 



and solving the RG equation for K d (d s ,fi R /^ 2 ,z\,Z2,^)- The constants K d (e) are identical 
toZ'^E) given in [43]. The constants gJ (e) are related to the finite boundary functions 
G d (a s (q 2 ), l,zi,Z2,e)- We define the Q d -(e) through the relation 



i=i 



^i' W (£) = E«^ 2 (l-Zi)(l-Z 2 ))^>) (30) 

i=l 



— / 
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The zi , Z2 independent constants Q d .(e) are obtained by demanding the finiteness of A s ^ 7 given 
in eqn.®. Without setting £ = in eqn.®, we expand A s ^ 7 as 

^dM^q 2 A^F^) = £4(A'I)A S J / l) (zi,Z2^ 2 ,^,^,e)- (31) 

Using the above expansion and eqn.© we determine these constants by comparing the pole as well 
as non-pole terms of the form factors, the mass factorisation kernels and the coefficient functions 
A^ 1 l ' expanded in powers of £ to the desired accuracy. Since the Gj(e) in the form factors are 

found to satisfy a specific structure in terms of / 7 ,|3 7 [43], we find that the constants Q di {z) also 
satisfy the following expansions containing these constants. 

?i«0 = -yf+E^Tjf . 

k=l 

OO 

V d ^) = -^-2Po?j[f + E^/f . 

k=l 

OO 

e ,>) = -A - - 2P» ffi; 2 (1) + 2Po? if + E ^i,f • 

/c=l 

?,(>) = -/i- 2 fc5j;, (1) -2Pl(ffi2' ) + 4Po?";.f ) ), 

Now that we have a better understanding [27] of the structure of even the single pole terms 
of the form factors, we can predict all the poles including the single pole of the soft distribution 
function from those of the form factors, the renormalisation constants and the mass factorisation 
kernels. The coefficients of the single poles are proportional to the constants —fj which are not 
only process independent but also maximally non-abelian. The £ dependent terms in Q d (e) can 
be obtained from the fixed order (in a s ) computations of cross sections and the finite parts of the 
form factors. At the moment, we know Q d l (e) to all orders in £, Q d 2 (e) to order £ and Q d 3 (e) 

to order e°. The lowest order term g dl (e) is known to all orders in £ from the exact fixed-order 
soft contribution at order a s . The next-to-leading order in a s (lowest order) computations of the 
total cross sections for DY and Higgs production determine the constants Q d j (e) and the results 
reveal that they are maximally non-abelian to all orders in £. One can similarly determine the £ 
dependent parts of soft cross sections beyond the order a s . The easier method is to use total cross 
sections that are known upto NNLO level. We can easily extract Q d 2 (e) upto order £ by using the 
fact that these constants are independent of Zj (j = 1 , 2) and the differential cross sections satisfy 
the relations 

jf 1 dx° x jf 1 dx\ [x\^f- X (4+4)^ = [ dx\f Q dx\ [x\x^f- X ^ = jf 1 dx x N - 1 a 1 , (33) 
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while the o are known for both DY and Higgs production upto NNLO level [17-27]. An alterna- 
tive method is to take /V — > °° on both sides of eqn. (l3*3l . In this limit, we find the following useful 
relation between the constants <j>^ (e) that appear in eqn.(l27l) and <j^'M(e) that contribute to the 
soft distribution function of the total cross section: 



^ (,> w=i^^'' (i) ( e >- < 34 > 



Both the methods give 



77/,(3) _ r f 1 r 2 



7,(1) Z2428 67 2 44 

+C I n f (- 3 ^ + ft; 2 + ^), (35) 



■ / 



where Q = Cf for I = q,b and C/ = Ca for / = g. Interestingly these constants Q d -(e) turn out to 
be maximally non-abelian. That is, they satisfy 

?i(e)=?£(e) = g?i(e). (36) 
This implies that the soft distribution functions for the differential cross sections satisfy 

<£/(a,,<,V,zi,z 2 ,e) =<Z> d b (a S: q 2 ,/j 2 ,zuz 2 ,£) = ^ 4>/ (a,,^ 2 ,^ 2 ,zi,z 2 ,£) , (37) 

wjcto orJer a 2 , similar to the soft distributions that appear in the total cross sections. We expect 
that this property will hold to all orders in perturbation theory because of the fact that it originates 
entirely from the soft part of the differential cross sections. 

The threshold corrections dominate when the partonic scaling variables z\ and zi approach 
their kinematic limit, which is unity, through the distributions 8(1 —Zj) and 'Di(zj) with j = 1,2. 
Resummations of threshold enhanced contributions are usually done in Mellin N space which has 
been a successful approach. See [28, 29, 60, 61] for resummation of total cross sections. We show 
in the following how the soft distribution functions 4>^(a J ,^ 2 , J u 2 ,zi,Z2,£) capture all the features 
of the iV space resummation approach. The exponents of the Zj (with j = 1,2) space resummed 
cross sections get contributions from the form factors through 8(1 — Zi)8(l — Z2) terms and from 
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the soft distribution functions through 8(1 — Zj) as well as the distributions 2>;(z/). We can rewrite 
the soft distribution function as 

*i(fi,,«V, Z i,za,e) = |s(l -«2) ^ T -l_.| j^ (1_ri) («,(X 2 )) 



+Gi(a,(^ 2 (l-z 1 )),£) 



+<; 



4(1-^(1-^1^ ~W 



+G; / (a,(^ 2 (l-z 1 )(l- Z2 )),£)}^) 



+18(1-^)8(1 -Z2) £4 (J) V e $/ W (e) 



4*>-*>(ry + £«@)V*c> 



where 



+(zi <->Z2)- 



G/(a,(^(zi,z 2 )),e) = £4 

(=i 



q 2 g(zi,Z2) 
M 2 



(38) 



(39) 



The third term in eqn.d3~8b contains the correct poles in £ to cancel those coming from the form 
factors as well as the 8(1 — Zj) parts of the mass factorisation kernels. The fourth term contains 
only poles that cancel against the D®(zj) parts of the mass factorisation kernels. The remaining 
finite terms as £ becomes zero in the first three terms contribute to the soft-plus-virtual parts of 
the differential cross sections. Hence, adding the eqn.OSb to the renormalised form factors and 
the mass factorisation kernels, performing the coupling constant renormalisation, and then finally 
taking the double Mellin moment mN\,N2, we get the resummed result analogous to the threshold 
resummation formula that one obtains for the total inclusive cross sections(see [28, 29, 60, 61]) 
when £ — > 0. A similar result for the resummed rapidity distribution scheme can be found in [29]. 

In the double Mellin space (N\,N2) the threshold enhanced differential cross section will be 
proportional to 



exp 



2 / dz\Z\ 




#1-1 



dz 2 z 2 - 9 dmts (a s ,q ,fi ,Zi,Z2) 



(40) 



Similar to the soft distribution functions Qp(a s ,q 2 ^fj 2 ,z,£) that enter in DY and Higgs produc- 
tion [42], the present <J>J (a s ,q 2 ,/u ,zi,Z2,£) are also maximally non-abelian. Using the resummed 
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result given in eqn.®, and the exponents g /(e) (see [35]), Q d ; (e), we can obtain the higher order 
soft-plus-virtual contributions to the differential cross sections. The available exponents are 

g{' J , Jdf for i' = a11 , 

82 J : jjf fOT 7=0,1 , 

■ ?£ W for j = , 

in addition to the known p\- (i = 0, 1,2,3), the constants in the splitting functions A^, B\ (i = 
1,2,3), the maximally non-abelian constants // (i = 1,2,3) and the anomalous dimensions "f?(i = 
0, 1 , 2, 3) . For I = q,g, the constants g|' 3 an ^ g\ 1 2X6 known for 7 = 2,3 also (see [34]). Using the 
resummed expression given in eqn.® and the known exponents, we present here the results for 
Aj V j for DY and Higgs production. Using our approach we have first reproduced the Drell-Yan 

coefficient A?' (l) , known upto NNLO (i = 0, 1,2) [50]. We then obtain A' v ; (l) and A^ (l) for the 

Higgs production upto NNLO (i = 1,2). For N 3 LO for I = q,b,g, a partial result A^ v } , i.e., a 
result without the 8(1 — zi)S(l — zi) part can be computed from our formula given in eqn.®. The 
coefficient of 8(1 — zi)8(l — zi) part depends on still unknown constants Cj\ , £3 • We 

can also obtain a result to N 4 LO order where we can predict partial soft-plus-virtual contributions 
containing everything except the terms in 2>o(z;)8(l — Zj), t Do(zi)£>o(zj),'Di(zi)&(l—Zj) and 8(1 — 
Zi)8(l — Z2) for the Drell-Yan N 4 LO coefficient ^„ , the gluon fusion to Higgs N 4 LO coefficient 

^dg 4 ^ anc * me bottom quark annihilation to Higgs boson N 4 LO coefficient A^ V £ . The results are 
presented in the Appendix B for = /jfi = q 2 '. The convolutions of distributions of the form 
£>l(zj) ® t D m (zj) for any arbitrary /, m can be done using the general formula given in [43]. Using 
these convolutions it is straightforward to calculate A*?} for i = 1, ...,4 for both DY (I = q) and 
Higgs (I = g, b) production. An alternative derivation of the NNLO DY soft-plus-virtual terms in 
the DIS renormalisation scheme can be found in [62], where they calculated di-lepton production 
cross sections at fixed target energies. 

The differential cross sections for I = q can be expanded in powers of the strong coupling 
constant as 

da 1 f ( . da 1 ^ 

~dY = Tjf'-dT- (41) 

We split the partonic cross section into hard and sv parts: 

do 1 ^ _ Jo hard ' / '(') do svJ 'W 
dY ~ dY + dY ' (42) 
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Figure 1 : Rapidity distributions for DY production at the LHC, and their (jlr scale dependence(with 
Mf = Q 2 = ^/ 2 +/-)- Here we denote M = Mi+ } -. The abbreviation "pSV" means partial-soft-plus- 
virtual. 



(43) 



The SM coefficients '® can be found in [50,51], and the 0%®$,%,$) for 

the Higgs are given in Appendix A. The soft-plus-virtual part of the DY partonic cross section can 
be expressed as 

j2 a sv,q,(i) r\ r\ 



x / dz\ j dz2 5(x? -xizi) 5(jc2 -X2Z2) Ay q ^\zuZ2,q l ^ F ^i) , 

'0 J 



(44) 
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Figure 2: Rapidity distributions for Higgs production through gluon fusion at LHC, and their /jr 
scale dependence for the Higgs boson of mass m H = 115 GeV with fi F = m H . 



and for the Higgs partonic cross section we have 

j a sv,g,(i) 



25 



dY 



Q H dx\ dx 2 ^ g (xi,x 2 ,i4) 



dzi dz2 5(jc? -xizi) 8(^2 -X2Z2) A^ g W {z\,zi,n? H ,i? F ,vfy . (45) 



X 



SV (l } / \ 

The coefficients A Y £ b \zi,Z2,q iMfiMr) are presented in the Appendix B, where we use the nor- 
malisation Ay^(zi,z 2 ,# 2 ,/4>/4) = 8(1 -zi)8(l -z 2 ). The constants g S M, q , Qh are given by 



QSM,i 



Aa 1 



Ql- 



2q 2 (q 2 -M 2 z ) 



* {^-M 2 Y + M 2 Y 2 )clsl 



n V V 
QqSeSq 



((q2-M 2 /+M 2 r 2 ) ci4 ( igVe)2 + ( ^ )2 ) ( ( ^ )2 + {gAq)2 ) 



nm 2 H G 2 B 
4(^V 2 -1) ' 



Qh = 



(46) 
(47) 
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The xf differential cross sections can be obtained from the Y differential cross sections by replacing 
D ab\AAil4) b y <i ( 'Vi>4>/4)- For the sv P^ we identify Ayf ] = (jt{ +JC§)~ 1 A^ i) with 
the replacement of D^^,^,^) b Y 

in the right-hand- side. The electro-weak 
constants appearing in eqns. (l46l47l) can be found in [25,51]. 

For our numerical results we choose the center-of-mass energy to be \/S =14 TeV for the LHC. 
The standard model parameters that enter our computation are the Fermi constant Gp = 4541.68 
pb, the Z boson mass M z = 91 . 1876 GeV and top quark mass m t = 173.4 GeV. The strong coupling 
constant 0C s (/4) is evolved using the 4-loop RG equations depending on the order in which the cross 
section is evaluated. We choose a\°(M z ) = 0.130, a^ LO (M z ) = 0.119, c^ NL0 (M z ) =0.115 and 
af LO (M z ) = 0. 1 14 for i> 2. We use MRST 2001 LO for leading order, MRST2001 NLO for NLO 
and MRST 2002 NNLO for N'LO with i > 1 [63, 64]. The impact of the soft-plus-virtual N 2 LO 
and the partial soft-plus-virtual N 3 LO contributions to the DY rapidity differential cross section 
at the LHC is presented in figure 1. Note that we have not plotted the partial soft-plus-virtual 
N 4 LO contributions because there are no N 3 LO parton densities. In the first plot we have shown 
the rapidity distribution in pb/GeV for a di-lepton mass of 1 15 GeV. For LO and NLO, we used 
the exact results which contain both the soft-plus-virtual as well as the regular hard contributions. 
For NLO (i = 2, 3), we use only the soft-plus-virtual results extracted from the resummed formula. 
Here we have set/j^ = /jr = 1 15 GeV. We find that the inclusion of NLO (i = 2, 3) terms only make 
small changes in the differential cross section which confirms the reliability of the perturbative 
approach. 

In the second plot of fig. 1 we have shown the scale variation of the rapidity distribution using 
the ratio: 

<tf = rt) _1 £(tf), (48) 

plotted as a function of /j//jq = /jR/\q\, where we have fixed fJ F = q 2 . It is clear from the second 
plot of fig. 1 that the inclusion of N'LO (i = 2, 3) soft-plus-virtual contributions further reduces the 
scale ambiguity. 

The impact of soft-plus-virtual parts N 2 LO and the partial soft-plus-virtual N 3 LO contributions 
to Higgs production through gluon fusion at the LHC is presented in figure 2. We see that the 
inclusion of the higher order terms does not make any appreciable change in the magnitude of the 
rapidity distribution. Again this confirms the reliability of the perturbation series. The second plot 
in figure 2 shows R in eqn.48 as a function of h/hq — /jr/wim, where we have fixed /j f = mjj, and 
demonstrates that the inclusion of the higher order terms reduces the sensitivity to the choice of 
the scale. 

We also present the numerical values of the rapidity distributions in figures 1 and 2 in Tables 
1 and 2 respectively. These numbers allow a more direct comparison with other theoretical papers 
and are useful to the experimental groups working at the LHC. 

Previous calculations of differential distributions in NLO using the effective Lagrangian (or 
the m t — > oo approach) have been presented in [65]- [69]. In the same approach the resummation 
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Y 


LO 


NLO 


N 2 LO S v 


N 3 LO p sv 


0. 


0.1326 


0.1735 


0.1734 


0.1732 


0.4 


0.1327 


0.1733 


0.1732 


0.1729 


0.8 


0.1330 


0.1729 


0.1725 


0.1722 


1.2 


0.1335 


0.1721 


0.1714 


0.1711 


1.6 


0.1343 


0.1707 


0.1698 


0.1695 


2. 


0.1346 


0.1678 


0.1670 


0.1667 


2.4 


0.1328 


0.1614 


0.1610 


0.1607 


2.8 


0.1273 


0.1483 


0.1485 


0.1482 


3.2 


0.1123 


0.1240 


0.1244 


0.1241 


3.6 


0.0832 


0.0858 


0.0857 


0.0855 



Table 1 : Values for d 2 c/dMdY (M=l 15 GeV) in pb/GeV at fixed values of Y which are plotted in 
figure 1. 



Y 


LO 


NLO 


N z LO S v 


N 3 LO p sv 


0. 


4.366 


7.236 


7.765 


7.848 


0.4 


4.294 


7.106 


7.632 


7.713 


0.8 


4.084 


6.722 


7.236 


7.313 


1.2 


3.731 


6.091 


6.585 


6.655 


1.6 


3.278 


5.264 


5.724 


5.784 


2. 


2.715 


4.265 


4.675 


4.724 


2.4 


2.068 


3.173 


3.512 


3.549 


2.8 


1.410 


2.087 


2.337 


2.361 


3.2 


0.779 


1.135 


1.285 


1.298 


3.6 


0.300 


0.448 


0.511 


0.516 



Table 2: Values for da/dY (for M H = 1 15 GeV) in pb/GeV at fixed values of Y which are plotted 
in figure 2. 



Y 


LO 


NLO 


N 2 LO S v 


N j LO p sv 


0. 


4.085 


6.791 


7.293 


7.370 


0.4 


4.016 


6.668 


7.166 


7.241 


0.8 


3.815 


6.298 


6.785 


6.856 


1.2 


3.477 


5.693 


6.161 


6.225 


1.6 


3.042 


4.899 


5.334 


5.389 


2. 


2.505 


3.944 


4.330 


4.375 


2.4 


1.888 


2.906 


3.223 


3.256 


2.8 


1.266 


1.883 


2.112 


2.133 


3.2 


0.680 


0.997 


1.131 


1.143 


3.6 


0.250 


0.376 


0.429 


0.434 



Table 3: Values for do/dY (for M H = 120 GeV) in pb/GeV at fixed values of Y. 
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of the logarithmically enhanced contributions have been carried out in [70] - [73]. Our results 
agree exactly with the DY NLO results for rapidity distributions in [19]. We cannot compare the 
NLO Higgs rapidity plots directly with those in [66] because there we used m# = 120 GeV and 
we had to impose a cut on the Higgs p t . However we have rerun our programs with m H = 120 
GeV to allow a better comparison with both the results in [66]. The number for the Higgs rapidity 
distribution are given in Table 3. Our numbers are also consistent with the normalized Higgs boson 
rapidity distribution in fig. 1 in [69]. These checks indicate that everything is consistent with the 
NLO results. 

We have compared our results for the Drell-Yan and Higgs rapidity distributions against the 
NNLO results published in [71,72]. Our soft-plus-virtual NNLO approximations agree very well 
with their exact NNLO results. Our partial soft-plus-virtual N 3 LO results are new and cannot be 
compared with any other calculation. 

To summarise, we have systematically studied the soft-plus-virtual corrections to differential 
cross sections in rapidity for DY and Higgs production through both gluon fusion and bottom 
quark annihilation. The resummation of these corrections has been achieved using renormalisation 
group invariance, Sudakov resummation of scattering amplitudes and the factorisation property of 
the hard scattering cross sections. Our analytical results are presented in Appendices A and B. 
It is now straightforward to obtain resummed threshold contributions to both xf and Y rapidity 
distributions of di-lepton pairs in the DY process and of Higgs bosons in Higgs productions. This 
requires a double Mellin transform in the space of two variables Ni and N2, see eqn.(40). Using 
our resummed results we have computed soft-plus-virtual differential cross sections at N 2 LO and 
partial-soft-plus-virtual differential cross sections in N 3 LO. Finally we have presented the numer- 
ical impact of these results on the rapidity differential cross sections. 

Acknowledgments: V. Ravindran and J. Smith would like to thank Prof. P. van Baal for 
hospitality. The work of J. Smith has been partially supported by the National Science Foundation, 
grant PHY-0098527. We all acknowledge support from the FOM under Grant No. L0102M and 
discussions with E. Laenen, S. Moch, J. Vermaseren and A. Vogt. 
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A Hart parts 



In this appendix, we list the C^'\x° l1 x ( 2 ) and D^' l \x c { 1 x ( 2 ) that contribute to the hard parts of the 
cross sections. We start by defining the following parton density combinations 

H q q(x h X 2 ^ 2 F ) = fq l {xi,l?p) ff(x 2 ^ 2 F ) + f£ (x h f? F ) f^ 2 (x 2 ,l4) , 

H gq (x h x 2 ,i4) = fg l {xuH 2 F ) (4 >2 (^2,A'f)+/|' 2 (^2,A'f)) > 
H qg (xi,x 2 ,n F ) = H gq (x 2 ,xi,iJ F ) , 
H gg {x u x 2 ,n 2 F ) = f^(x h /4) fg 2 (x 2 ^ 2 F ). 



(1) 



In terms of these combinations, we list the (AA) that 

appear in the hard parts of the x F - 

differential cross sections 

Cgg^ \AA) = HggiAiAiA) i 

H gg (xi,X 2 ,fi F ) 



Cg 8 (l \AA) - C \j^ "4+1) 



,xf_ ,A_%_ 4\ 

3 ~^ 2 ~^J3 )^ a \ 

X\ XTy J 



4 



{xx- A) 



L 



, 1 Hssiixux^.ul) 
+ n dxi 88 ,u ' 2 ^ F) 



•1 r\ 

+ I dx\ / „ dx2 



x\ ( Xl - 

Hg gA (x U X2,^ 2 F ) 



J x^ ) ("^1 - ^ - '^) 



— 4 -j- + 4 -| 1 



+ 



X] * I4 



dX2 



2H gg ^ 2 {x h X2,n 2 F ) 
( Xl -^)( X2 -x\){x\+x%) 



Hgg(x\,X2,Hp) 



2L (jfj , x\) log + (x $ ,^)) 2 + 6^ 2 



+ 



Hgg(xi,x 2 ,Li F ) 



-~ ( Ax\x\ x\ 
xi \ 



l£+A) 

Sx\ lx% (XI+X2) 3 (XI+X 2 )(X2+X 1 )(X 1 +X° 2 ) 

+4x 2 2 (xfxl +x° 2 xf) +4x 2 (A 3 A 2 +A 2 A 3 )) + ^ (^xlixf -^x\) 

+x 2 2 (4xf + 12*° 2 A -*AA 2 ) +X2(l6xfx° 2 + Uxfxf -A^xf) 
+&{xfxf+xfxf)\ + — (x 3 2 (Sx° 2 - Ax\) +xl(Sxf + Sxf - Sx ^) 
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+x 2 (8x° 2 3 + 12xf + 8^ V 2 - foP^ 2 ) + 2^(12x5 V 2 2 + 12xfxf ) 
+ 12^V 2 + 20^V 2 2 + 8^ 2 3 ^ +*?(l2-8^ 



+x\ (l4x\ + 4^- - 4x° 2 ) +xi (^x° 2 + 36xf - 4^- +x 2 ( - 4i 
+20x* - 4- + 20x° 2 ) +4(-A-l2^ + 4o) + 4 (4 ~ 4 



+ 10x° 1 r° 2 2 + lOxf x^ + 2xf + 2xf 



+ (1-2), 



Cf{ k dxi (A+A) 



2^j 4 4 



X| a*i .xfj y 



+ 



+ 



1 ^ H qg , 2 (x h X2,ti) 



ip^l — -*"2) +^2) 

/n L dx 2 H qg (xi,x 2 ,i4) 



2w 2^ 4 4' 



\ x x x] 



1 4 4 

(^+x° 2 )(x 2 +^) V X^ XI x\ 



2 4 
+ — o + 



X 2 X2 xix 2 (xi+X2) 2 \ ijxj x\x\ x 2 



X\ x 2 



1 



(xi +X2) V *2 



1++2 



^ dxi L 

A^ J ^2 



(A +A) 



1 



(xi +x 2 ) 3 \xj 

JCs) n A^j 



2 (4x^ + 4^j Xj) 



+4^-1^+4^-1^) + 



J! 



1 Y4+ 1 



(xi +x 2 ) Vx°i -4 



+ (1-2). (2) 
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We have introduced the following abbreviations 

H a b,12( X hX2,&) = H ab {xi,X2,l? F )-H ab {^X2,l? F )-H ab {x\^ 2 ,l?p) 

+H ab (x^,x^,/d F ) , 

H ab:1 (x h Z,IJ F ) = Ha^Xi^^-Hah^Z^p), 

H ah ^z,X2,i? F ) = H ab (z,X2,M F )-H ab (z,x° 2 ,iii F ), 



In 



In 



q 2 (A+4)(i-4)(xi-A) 

1+2 ^ £(x^) = ln 



L bl = In 



q i(l-xl)( Xl -x\) 



(l-x*)(l-xl) 



The D^ l \x Q l ,x° 2 ) that appear in the hard parts of the rapidity distributions are listed below 



H gg;1 (x h x° 2 ,n 2 F ) 



2 .v, ' ,'|J' A '"- 



J^j (-^1 -^j ) 



4%, 



+ 



+ 



- 1 H ggA (x h X 2 ^ 2 F ) 

1x9, ( X1 -x\) 



1 dx 2 2H SsM X ^ X 2^f) 



■ 4— + 4— 1 

XT tP 



li dXi L 

lx\ ^ lx% ^ ( Xl -^)(X2-X\) 

+H gg (x° v x° 2 ^ 2 ) 2*(j°,a°)log^ +6C2+(3CW4)) 



<&2 



H gg (xi,x 2 ,i4) 



1 ^JJ7 



4-x 2 x^ ~\~ Sx 2 x^ $ 2 



+ ix] ^ lx\ (xi +X° l )(x2+d 2 )(xiX 2 +X l X 2 ) 4 

+8x 2 2 x° l 7 xf + Sx 2 xfxf^ + 1 V 2 4 + 16x1^^ + 32^x° 1 6 x° 2 
+24x 2 x 1 6 x° 2 3 ^ + 1 (^Oxfrf + Ax\x\ 5 + Sxl^xl + 32x 2 ^ 5 xf 



(3) 



(4) 
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+4 wj 5 xf + y/"^) +4 ( 4 |r) +4 ( 4 ^ 4 ) +4 (i24A4 



+2Sx° l xf }+x 2 l ( 48r] 2 ^ 4 + I6xixf + 2Qx\xfxf + 88*2*?^ 



\2 n3 



+ X1 (^Sxfxf + + Ax\^ + Sxlxfxl + 64x 2 xfxf 



x\ 



+40a° 1 > 2 4 



+ (1-2), 



CF |j^o dx\H qg {x\,x\,i? F ) 



2r} 4 4 \ 2x? 

Xj^ y 



H qga (x h x 2 ^ 2 F ) 



(X2-X° 2 ) 



2 ^ r 2x°, 4 4 
1 + 



X 



2 XI A? 



■1 rl 

+ I dx\ I dx2 



H qg (xi,X2,fIf) 



'jfj Jx\ ~* ( Xl +X° 1 )(X2+4)(X 1 X 2 +X° 1 X2^ [x\\ 21 
— 4-X 2 X^ X^2 — 4-X2X^ X^ ^ H ^^.JC^^j — ^^2^\ 4. — 4 ^2^j 4. ^ 

6x2X^ X^2 ~\~ 2xj X^2 — 12x2X^X^2 ~\~ 4^2-^j — 2XjX^ J 



^/ \4^4) 



l<->2 



Cj- <( / n Jxi / n dx2 



H qq (x h x 2 ^ 2 F ) 



l/„ 6n3 



V] (x 1 +X Q l )(x 2 +X 2 )(x 1 X Q 2 +X l X 2 ) 4 

+Sx° l 6 x° 2 4 ) + - (l6x 2 xfxf + Sxl^xf + Sxfxf) +jci ( - 32x 2 xfx° 2 



Xl 



21 



+ (1-2), 

where we have introduced the following abbreviations 

V(1-jS)(*1-jPi) 



%a x = In 



=ln 



q 2 (l-A)(xi-A) 



In 



2i] 



The %a 2 , and 3Cc 2 can be obtained from %a x , Kj >l and «Ccj by using 1^2 symmetry. 



(5) 



(6) 



B Soft-plus-virtual parts 

Here we list below the Ay V ^ that contribute to the soft-plus -virtual parts of the cross sections for 
the choice f? R = fip = q 2 . 



v sv,(l) _ 



A sv,(2) 



S(l-zi)8(l-z 2 ) 

+ (Z\ — Z2), 

8(l-zi)8(l-z 2 ) 



Cf 



-8 + 6£ 



+ © ©o 



Cf 2 



127 76 4 
«/Cf ( - - ~g & + 3 & ) + C F C A 



+ 2>i8(l-zi) 



1535 



C F 4 



24 



430 26 2 86 \ 2 /511 152 2 

"~9~ ^-yC 2 + yC3)+C F [— -67^2+^-^2-30 C 3 
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+ £> £>1 

+© 8(l-zi) 
+C^32^ 3 
+0^-64+16^ 
+© 3 8(l-zi)[c|| 8 



n f c n -^) +CfCa \ 9 
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iifCr[ - I + C f Ca( -y 



4^2 +Cj -32 + 8 ^ 



+ 2>i2)i 



+ £>0©2 



Cj| 24 



C|| 24 



nfCr[ lll.^ 2 ) + c F c A (-^ + 4 4k + 2H, 



27 



+ ©i8(l-zi) 



27 3 
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+ © 2 8(1-zi) 



"8^2 

»/C>( 3 J +C F C A I -y 
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+ {Z\ <->Z2), 



2) 2)o 



4102 256 ^ , 

n f C F C A - + — ,2 J + n f C F 



81 



536 224 „ 160 „ 
H ^C2 + -^^3 



+«fC F 



200 16 
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/15503 _340 88 

81 3 s 5 S2 Si 
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+C F C A 
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/2312 32 
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+ ©0©2 

/1072 



h/CfCa 
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